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ABSTRACT 


Estimation  procedures  for  the  shape  and  scale  parameters  of  a 
mixed  exponential,  decreasing  failure  rate  distribution,  are  compared 
and  the  more  advantageous  methods  are  determined  for  various  situations. 

When  the  scale  parameter  is  known  and  the  shape  parameter  is  to 
be  estimated,  a  transformation  of  the  data  reduces  the  problem  to 
that  of  estimating  the  failure  rate  of  an  exponential  distribution, 
in  which  case  the  Maximum  Likelihood  Estimate,  MLE,  involves  the 
total  life  statistic,  the  properties  of  which  are  well  known. 

A  Best  Linear  Unbiased  Estimator  of  the  scale  parameter,  based 
on  the  first  k  order  statistics,  BLUE-k,  is  more  accurate  than  the 
MLE  in  the  case  of  type  I  or  type  II  censorial  for  small  sample  sizes 
but  is  shown  to  be  generally  less  accurate  for  random  censoring  of 
a  type  often  found  in  industrial  applications. 

In  the  case  both  p« raneters  are  laUmown,  ;n<  BUS  methods  are 
not  available.  As  is  usually  the  case,  the  joint  lit,  despite 
being  only  implicitly  defined,  are  much  superior  in  performance  to 
estimates  made  using  the  method  of  moments.  In  fact,  the  distribu¬ 
tion  function  estimated  using  the  joint  MLE's  of  the  parameters  is 
surprisingly  closer  to  the  true  distribution  for  regions  of  interest 
in  reliability  theory,  than  is  the  estimated  distribution  function 
using  BLUE-k  for  the  scale  parameter  and  a  known  shape  parameter. 


1.  INTRODUCTION 


A  mixed  exponential  distribution,  sometimes  called  the  Lomax 
or  Pareto  distribution  of  the  second  kind  has  for  parameters 
o, 0  >  0,  the  representation 

F(x)  -  1  -  (H-x/8)"*  for  x  >  0.  (1) 

This  distribution  has  proved  to  be  useful  for  representing  the  failure 
times  of  electronic  packages ,  Myhre  and  Saunders  [1979] . 

The  estimation  of  the  parameters  of  this  distribution  has  been 
discussed  in  numerous  papers.  Kulldorff  and  VSnnman  [1973]  have 
derived  the  best  linear  unbiased  estimate  (BLUE)  of  the  scale  para¬ 
meter,  0,  for  known  values  of  the  shape  parameter,  o  >  2.  For 
known  a  <  2  an  asymptotic  best  linear  unbiased  estimator  (ABLUE) 
for  8  was  also  derived.  Both  of  these  procedures  require  a 
complete  sample.  Vannman  [1976]  has  obtained,  for  known  a,  a  BLUE 
estimate  of  8  from  censored  data  (BLUE-k) .  The  only  restriction 
on  the  procedure  is  that  k<n+l-2/o  where  n  is  the  sample 
size  and  k  is  the  number  of  failures.  Each  of  these  estimation 
procedures  assumes  a  knowledge  of  the  shape  parameter,  a  >  0. 

Whenever  the  scale  parameter  8  is  known  the  transformation 
ln(l+x/8)  makes  the  life  observed  exponentially  distributed  with 
unknown  failure  rate  a,  in  which  case  the  optimum  property  of  the 
total  life  statistic  for  type  I  or  type  II  censoring  is  well  known. 
Consequently  we  do  not  discuss  the  case  further. 
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Harris  and  Singpurwalla  [1968]  suggested  the  method  of  moments 
as  an  estimation  procedure  wherever  S  is,  or  a  and  8  are,  un¬ 
known.  This  method,  however,  requires  not  only  a  complete  sample 
but  the  restriction  a  >  2,  for  otherwise  the  moments  fail  to  exist. 
The  only  point  in  favor  of  this  method  of  moments  seems  to  be  that 
it  gives  simple  closed  expressions  for  the  estimates  of  both  a  and 
6.  Unfortunately,  we  find  that  they  fail  to  exist  so  often  in 
practice  as  to  be  virtually  unusable.  For  examqle,  in  our  simula¬ 
tions  using  complete  samples  as  large  as  n  ■  25  with  a  *  4,  the 
estimates  fail  to  exist  nearly  30%  of  the  time. 

Harris  and  Singpurwalla  [1969]  first  exhibited  the  maximum 
likelihood  equations  of  this  family  for  complete  samples.  In  Myhre 
and  Saunders  [1979]  maximum  likelihood  equations  are  derived  for 
both  the  shape  parameter  o,  and  the  scale  parameter  g,,  for 
censored  or  complete  samples. 

While  the  maximum  likelihood  estimates  are  only  implicitly 
defined  as  the  solutions  to  certain  equations,  this  causes  little 
practical  difficulty  since  they  can  be  calculated  with  a  programmable 
hand  calculator. 

For  given  a,  the  MLE  for  6,  (say  6),  is  shown  by  Myhre 
and  Saunders  [1979]  to  be  the  solution  of  the  equation  A(l/B)  ■  0, 
where 

n  .  .  k  . 

A(x)  «  o  l  (1+t.x)  ♦  l  (1+t.x)"1  -  no  , 

j«l  3  1  3 
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with  denoting  the  k  failure  times  and  t^,...,^ 

denoting  the  n-k  success  times  in  the  sample  (tlf...,t  ). 

Since  A(o)  ■  k,  A(«)  »  -no  and  A(x)  is  a  decreasing  function 
in  x,  a  unique  solution  always  exists  for  any  positive  o  or  any 
type  of  censoring.  This  is  in  contrast  with  the  BLUE  which  requires 
a  complete  sample  and  a  >  2,  the  ABLUE  which  needs  a  complete 
sample  and  the  BLUE-k  for  which  k  <  n+1  -  2/a  must  hold.  (Note 
that  for  modest  n  and  small  a,  k  is  very  restricted.) 

In  this  paper  the  MLE  procedure  of  Myhre  and  Saunders  [1979]  is 
compared  with  the  BLUE  and  ABLUE  procedures  described  in  Kulldorff 
and  Vannman  [1973] .  For  complete  samples  with  a  >  2  we  know  that 
the  BLUE  and  the  MLE  have  an  asymptotic  relative  efficiency  of  one. 
However,  in  reliability  applications,  the  sample  is  often  censored 
and  some  data  from  industry  indicates  that  a  <  2  can  arise.  For 
these  cases  the  variance  of  the  MLE  Has  not  been  represented  in  a 
closed  form  and  hence  no  relative  efficiency  with  the  BLUE-k  can  be 
computed.  But,  the  asymptotic  variance  of  the  MLE  is  known  in 
these  cases,  as  was  shown  by  Myhre  and  Saunders,  [1980],  In  the 
appendix  it  is  shown  that  for  type  II  censoring  with  large  n  and 
k,  the  variance  of  the  MLE  is  less  than  the  variance  of  the  BLUE-k. 

For  small  n  and  k,  realistic  comparisons  are  made  using  simulation 
techniques  with  particular  attention  given  to  censored  data  for  the 
case  with  a  S  2. 
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2.  SIMULATION  MODEL 

To  compare  the  BLUE  and  the  MLE  of  8  for  fixed  a,  ordered 
failure  times  are  generated  from  the  distribution  in  equation  (1) 
with  6  ■  1.0  and  the  two  estimates  are  computed.  These  simulations 
were  run  for  n  *  25,  50,  100,  200;  a  ■  2.2,  2.4,  2.6,  2.8,  3.0, 

4.0.  For  a  fixed  o  and  one  sample  size,  one  hundred  simulations 
were  run  and  from  these  date,  the  averages  of  each  of  the 
two  estimates  of  8  and  their  sample  standard  deviations  about 
0  ■  1.0  were  computed.  It  was  seen  that  the  comparisons  did  not 
change  significantly  if  more  simulations  were  run  for  each  case,  so 
100  was  used  to  keep  the  computer  time  low.  For  complete  samples, 
with  a  >  2,  the  BLUE  was  seen.  Table  1,  to  be  no  better  than  the 
MLE. 

The  ABLUE  introduced  for  the  case  o  <  2,  is  a  more  difficult 
procedure  to  use  than  the  BLUE.  In  simulations  comparing  the  ABLUE 
and  the  MLE,  ordered  observations  are  generated  for  a  known  a,  and 
both  ABLUE  and  MLE  estimates  of  8  are  computed  from  this  complete 
sample  and  compared.  The  ABLUE  is  based  on  ten  order  statistics, 
the  largest  number  for  which  computational  results  are  available 
in  Kulldorff  and  Vannman  [1973] .  These  simulations  are  run  for 
n  »  25,  50,  100,  200;  o  ■  .5,  1.0,  1.5,  2.0,  3.0,  4.0.  Again 
for  a  known  a,  and  one  sample  size,  one  hundred  simulations  were 
run  and  from  these  simulations  the  averages  of  each  of  the  two 
estimates  of  8  and  their  sample  standard  deviations  about  0  ■  1.0 
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were  computed.  The  MLE  was  a  closer  estimator  of  g  •  1.0  as  is 
seen  in  Table  2. . 

The  estimation  procedures  developed  in  Vannman  [1976], 
(referred  to  as  BLUE-k) ,  produce  a  more  versatile  estimator  than 
the  BLUE  and  ABLUE  developed  earlier  (loc.  cit.) ,  and  therefore  we 
subjected  it  to  a  mere  thorough  investigation.  In  order  to  reflect 

practical  applications  in  which  tests  are  not  only  truncated  on 
time  (Type  I  censoring]  or  number  of  failures  (Type  II  censoring), 
but  are  also  randomly  truncated,  simulations  were  conducted  for 
all  these  types  of  censoring.  The  BLUE-k  is  based  on  k  order 
statistics,  where  k  is  known.  It  is  possible,  however,  to 
compute  an  estimate  of  g  based  on  the  BLUE-k  procedure  using 
randomly  censored  data  by  letting  k  be  the  observed  number  of 
failures;  while  this  procedure  is  not  theoretically  justified,  it 
is  a  natural  way  to  attempt  to  use  the  estimator  in  such  cases. 

The  simulation  comparing  the  BLUE-k  and  the  MLE  of  g  was 
again  run  for  Type  I,  Type  II  and  random  censoring  with  four  sample 
sizes,  n  *  25,  50,  100  and  200.  Samples  that  illustrate  the 
results  for  Types  I  and  II  censoring  are  grouped  in  Tables  3  and  4. 
For  these  samples,  a  runs  through  a  selected  set  of  values  from 
0.1  to  4.0.  The  pattern  that  emerges  does  not  depend  on  the  type 
of  censoring:  the  averages  and  sample  standard  deviations  for  all 
methods  improve  with  sample  size,  especially  for  the  MLE;  the 
BLUE-k  is  best  for  small  a,  a  ■  .1,  for  all  sample  sizes  checked 
but  the  difference  is  small  when  n  ■  200.  For  slightly  larger  a. 


a  •  .2,  .2,  .4,  the  BLUE-k  does  better  than  the  MLE,  i.e.,  has 
less  bias  and  smaller  sample  standard  deviation,  when  n  <  SO,  but 
the  results  are  very  close  when  n  *  100  and  n  «  200.  For 
larger  o  the  results  for  n  »  100,  200  are  very  close  with  a 
slight  edge  to  the  BLUE-k. 

Often  in  industrial  applications,  on  doesn't  find  Type  I  or 
Type  II  censoring  but  instead  encounters  random  censoring.  Because 
of  this,  simulations  with  random  truncation  times  were  run.  The 


censoring  distribution  function 


Cy+d-t)2 

(m+l)u 

l-(T*mD-v): 

mCm+DD* 


QCy) 


when  T-D  <  Y  <  T 


when  T  «  y  <  T+mD 


for  positive  m,  D  and  T,  was  chosen  for  its  versatility  and 
simplicity.  Here  m  >  1  makes  the  distribution  skewed  to  the  right 
and  m  <  1  makes  it  skewed  to  the  left.  The  spread  of  the  truncated 
times  is  given  by  (m+l)D.  Using  this  censoring  distribution  for  a 
fixed  T,  0,  and  m,  a  truncation  time  is  generated  for  each  item 
an  test.  If  the  generated  failure  (observation)  time,  t^,  is 
greater  than  y^»  then  the  item  is  a  success  and  t^  is  set  to 
Y^,  while  if  t^  2  it  is  counted  as  a  failure  and  is  unchanged. 

The  simulations  using  random  censoring,  compared  the  BLUE-k 
with  the  MLE  of  8,  and  were  run  with  sample  sizes  of  25,  50,  100 
and  200.  For  each  specifie .  o,  different  values  for  D,  m  and 
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T  were  specified.  D  ranged  from  .ST  to  0.8T,  m  >  1,2,3,  and  T 
from  F"1(0.S)  to  F"*(0.8)  where  F  is  given  in  equation  (1). 
These  simulations  showed  the  MLE  to  be  better  than  the  BLUE-k  as 
an  estimator  of  6  except  for  small  a  coupled  with  small  sample 
sizes.  (For  a  <  1.0,  n  ■  25  the  BLUE-k  has  less  bias  than  the 
MLE.)  These  are  the  cases  where  the  random  censoring  had  little 

chance  to  show  its  effect.  But  for  SO  5  n  5  200  and  o  >  0.2 
the  MLE  showed  itself  to  be  more  accurate.  These  results  follow 


in  Table  5. 


3.  JOINT  ESTIMATES 


The  comparisons  so  far  have  been  made  for  estimates  of  B  when 
a  is  known.  For  many  applications  in  reliability  the  parameter  a 
must  also  be  estimated.  The  joint  MLE  of  o  and  B  for  a  sample 
consisting  of  k  failure  times  t^,...,!^,  and  n-k  success  times 
tls+l,****tn  re<*u*res  a  numerical  solution  for  x  in  H(x)  •  0 
where 

H(x)  *  1/k  l  (1+t.x)-1  l  InCl+t.x)  -  l  r— . 
i-1  1  j-1  J  j-1  1  'jx 

n 

We  obtain  $  »  1/x  while  a  a  1/k  £  ln(l+t.x) .  A  necessary  and 

i*l  1 

sufficient  condition  that  the  solution  exists,  viz.,  H(x)  *  0  for 
some  x  >  0,  is  that 


If  this  condition  does  not  hold,  the  distribution,  F,  is  assumed 

»  n 

to  be  exponential  with  estimated  failure  rate,  \  ■  (1A)  1  t-. 

i*l  1 

No  BLUE  methods  exist  to  estimate  a  and  6  jointly.  Harris 
and  Singpuwalla  [1968]  used  the  method  of  moments  assuming  a  >  2 
as  an  estimation  procedure  for  a  and/or  8  but  only  for  complete 
samples. 

Ihe  joint  MLE  will  first  be  compared  with  the  BLUE-k  with 
a  assumed  known.  Then  the  joint  MLE  is  compared  with  the  joint 
estimation  by  method  of  moments. 

These  comparisons  will  be  made  using  the  distance  of  the 
estimated  distribution  from  the  true  distribution  in  the  region  of 


9 


interest,  not  by  using  a  measure  of  the  estimated  a,  8  from  the 

true.  This  is  thought  to  be  the  more  natural  method.  Each 

«• 

estimated  distribution  F  is  compared  to  the  true,  F,  at  ten 
evenly  spaced  values  from  the  interval  C0,F*1(.2))  by  the  stmt 

10 

S  »  I  |FCx.)  -  F(x.) [  . 
i-1  1  1 

This  is  computationally  easier  than  finding  the  maximum  distances 
between  the  functions  over  the  same  interval.  Tables  6  and  7  contain 
the  results  comparing  the  estimated  distribution  F^  using  MLE 

m 

estimates  of  o  and  8,  with  the  estimated  distribution  Fg 
with  the  true  a  and  the  BLUE-k  estimate  of  8*  When  the  data 
dictates,  F^  is  an  estimated  exponential  distribution. 

From  Table  6  (Type  I  censoring)  it  is  seen  that  for  the  cases 
examined,  the  accuracy  of  the  estimate  of  the  Mixed  Exponential 
distribution  function  using  joint  MLE's  for  unknown  a  and  8  is 
surprisingly  better  than  the  accuracy  of  the  estimated  distribution 
obtained  with  a  fixed  and  using  the  BLUE-k  estimate  of  8.  Under 
Type  II  censoring  the  accuracy  of  the  two  estimation  procedures  appears 
to  be  approximately  the  same  even  though  a  is  assumed  known  when 
using  the  BLUE-k  estimate  of  8. 

To  get  a  fairer  comparison  of  the  behavior  of  the  joint  esti¬ 
mates  of  a,  8,  the  parameter  a  must  also  be  estimated,  and  the 
only  alternative  procedure  is  the  method  of  moments  which  restricts 
one  to  complete  samples  with  a  >  2.  Even  for  complete  samples 
of  size,  say  n  »  25  or  SO,  the  method  of  moment  estimator  fails  to 


be  admissable  a  significant  ntmber  of  times.  In  Table  8,  some 
samples  are  given  comparing  the  joint  MLE  with  joint  method  of 
moments  estimates. 

Thus  for  a  very  important  case  in  applications,  namely 
censored  data  with  both  parameters  unknown,  only  the  maximum  like¬ 
lihood  estimates  are  available  and  they  are  seen  to  provide  reasonably 
accurate  estimators  of  the  distribution  in  the  lower  percentiles, 
the  region  of  our  interest. 
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4.  CONCLUSION 

For  the  case  where  the  shape  parameter  o  >  2  is  known,  with 
complete  samples,  the  BLUE  of  the  scale  parameter  0,  was  no  better 
than  the  MLE.  For  known  a  >  2,  with  complete  samples  the  comparison 
between  the  ABLUE  and  the  MLE  showed  the  MLE  to  be  the  better 
estimator.  For  censored  samples,  moreover  simulations  show  that  the 
MLE  is  as  accurate  an  estimate  of  the  scale  parameter  0  as  the  BLUE-k 
except  for  small  sample  sizes  and  for  small  values  of  the  shape  para¬ 
meter,  a. 

In  practice,  both  parameters  must  be  estimated  jointly.  Such 
joint  estimation  is  not  possible  with  BLUE,  ABLUE,  or  BLUE-k  methods, 
and  the  method  of  moments  fails  too  often.  Therefore,  in  applications 
the  joint  MLE  procedure  would  be  used  even  if  the  MLE  of  B  were  not 
more  accurate.  Moreover,  the  estimate  of  the  distribution  function 
using  joint  MLE's  of  the  scale  and  shape  parameters  compares  favorably 
or  looks  superior,  in  the  regions  of  interest  in  reliability  theory, 
to  the  estimate  of  the  distribution  function  obtained  using  the  BLUE-k 
for  the  scale  parameter  and  a  known  shape  parameter. 
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APPENDIX 

Let  e  >  0  be  given  and  denote  g  and  g  as  the  BLUE-k  and  the 
MLE,  respectively,  of  g.  Then  from  Vannman  [1976],  we  have 

VM!^  ■  ■  m 

where  R  is  the  ratio  of  gamma  functions,  namely 


R  -  r(n+l-<i») r(n-k+l)  /r(n+l-k-w)  r(n+l) 


with  u  »  2/a. 

From  Halperin  [1954],  the  asymptotic  variance  of  g  is  given  by 


AVAR(g) 


g2(q+2) 

iwtl-Cl- 


from  n  large. 

If  k  *  qn  for  some  q  e  (0,1)  and  n  -  k  >  u,  we  will  show  that 
the  ratio  of  variances  satisfies 


AVARfl)  1  '  “n  ~  £>»  .  t 

VAR(S)  1-(1-<0I4“ 

with  equality  of  u  »  0  or  in  the  limit  as  n  ♦  <*. 

Using  Sterling's  approximation,  we  find  for  the  ratio  R 

T  «/(■<-«)  (n-10(n-.),,-M(..k)n-k 
»  (n-w-k)n  (n-u-k)n”“  nn 

The  square  root  obviously  approaches  unity  as  h  -*■  «,  while  the  second 
can  be  rewritten 
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which  can  be  seen  to  approach  (l-q)w  since  k  »  qn.  Thus*  in  the 
limit  as  n  *►  •,  we  see  quality  is  attained. 

Me  want  to  exaaine  the  derivative  of  the  ratio  of  variances  as  a 
function  of  n  >  «  +  k.  Since  the  denominator  is  independent  of  n  we 
find  the  derivative,  call  it  D,  to  be  proportional  to 

D  -  X-  (1-R)  -  (1-q-  £)R(ln  R)  ’ . 
n2  n 

Since  0  <  R  <  1,  it  is  sufficient  to  show  only  that  (In  R)  •  <  0,  to 
conclude  D  >  0  from  which  it  follows  that  the  ratio  of  the  variances 
increases  to  unity.  We  check  that,  setting  p  -  1-q, 

(In  R)'  ■  iKn+l-w)  -  tKn+1)  ♦  pfrKnp+l-w)  -  ♦(np+1)] 

where  the  psi  (digamma)  function  is  the  derivative  of  the  logarithm  of 
the  gamma  function,  see  Abramov it z  and  Stegun  [1964] .  And  since  the 
psi  function  is  increasing  for  non -negative  argument,  loc.  cit.,  the 
result  is  proved. 
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Table  1 


IS 


Simulation  experiments  (100  realizations)  coopering  MLE  and 
BLUE  estimates  of  scale  parameter  6  *  1.0. 


Average  of  Estimates 


n 

«  25 

n 

N 

in 

© 

n  « 

100 

n  ■ 

200 

a 

MLE 

BLUE 

MLE 

BLUE 

MLE 

BLUE 

MLE 

BLUE 

2.2 

1.02 

1.00 

1.01 

1.00 

1.01 

1.00 

1.00 

1.00 

2.4 

1.03 

1.00 

1.04 

1.02 

1.01 

1.00 

1.00 

.99 

2.6 

1.07 

1.05 

.99 

.98 

1.01 

1.00 

1.00 

1.00 

2.8 

1.02 

1.00 

1.02 

1.00 

.99 

.99 

1.00 

1.00 

3.0 

1.02 

1.00 

1.01 

.99 

1.00 

.99 

1.01 

1.01 

4.0 

1.06 

1.04 

.99 

.98 

.99 

.99 

.99 

.99 

Sample  Standard  Deviations 


25 


50 


100 


MLE  BLUE 


MLE  BLUE 


MLE  BLUE 


2.2 


.28 


.27 


20 


.19 


.15 


.15 


2.4 


.33 


.32 


19 


.18 


.16 


.16 


2.6 


.33 


.32 


18 


.18 


.13 


.13 


2.8 


3.0 


4.0 


.27 


.26 


.26 


.27 


.12 


.12 


.25 


13 


.13 


.26 


.12 


.12 


n  »  200 
MLE  BLUE 


.09  .09 
.10  .09 
.10  .10 
.10  .10 
.10  .10 
.08  .08 


17 
16 

18 


16 

16 

18 


16 


Table  3.  Simulation  experiments  (100  realizations)  comparing  MLE  and 
BLUE«k  estimates  of  the  scale  parameter  8  ■  1.0  under  Type 
I  censoring  of  data. 


Averaae  of  Estimates 


BLUE-k  MLE  BLUE-k 


1.23  .95 


1.22  1 


1.17  .98 


1.21  1.02 


1.49  1.02 


1.32  .92 


-  100 


BLUE-k 


1.08  .91 


1.12  .99 


.98 


1.02  .95 


n  »  200 


MLE  BLUE-k 


.92 


1.03  .98 


1.04  1. 


1.01  .97 


1.06  1.02 


.87 

1.10  1. 

00 

1.04  1. 

1.18  1.01 


1.20  1.02 


1.17  1.02 


1.11  .98 


1.10  1.01 


1.08  1.02 


1.08  1.01 


1.10  1.02 


.98 


1.0S  1.01 


1.03  1. 


1.04  1.01 


0 

B9K9EH 

2  .99 

1.02  1. 

1.35  1.04 


1.21  .92 


1.18  1.05 


1.06  .94 


1.05  .98 


1.10  1.03 


le  Standard  Deviations 


n 

.  25 

n  ■ 

■  50 

n 

-  100 

n 

»  200 

MLE 

BLUE-k 

MLE 

BLUE-k 

MLE 

BLUE-k 

MLE 

BLUE -l 

3.67 

1.63 

1.04 

.75 

.49 

.41 

.31 

.39 

.64 

.80 

.51 

.35 

.28 

.26 

.25 

.49 

.71 

.50 

.39 

.33 

.20 

.19 

.88 

.65 

.60 

.49 

.35 

.32 

.21 

.20 

1.35 

.79 

.65 

.48 

.29 

.28 

.21 

.20 

1.12 

.57 

.38 

.35 

.42 

.27 

.20 

.18 

.92 

.63 

.56 

.40 

.31 

.28 

.18 

.17 

1.23 

.67 

.62 

.38 

.32 

.28 

.18 

.17 

1.16 

.68 

.51 

.44 

.30 

.26 

.19 

.19 

1.42 

.67 

.58 

.44 

.31 

.26 

.18 

.17 

.71 

.44 

.40 

.32 

.26 

.24 

.18 

.17 

KKH 

.53  .41 

.29 

.26 

.19  .19 

.77 

.43 

.48 

.37 

.34 

.25 

.17 

.16 

18 
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Table  4.  Simulation  experiments  (100  realizations)  comparing  MLE  and 
BLUE-k  estimates  of  the  scale  parameter  8  ■  1.0  under  Type 
II  censoring  with  k  ■  ,2n. 


Average  of  Estimates 


Table  5.  Simulation  experiments  (100  realizations)  comparing  MLE 

and  BLUE-k  estimates  of  the  scale  parameter  8»  1.0  under 
random  censoring. 


n  ■ 

>  25 

n 

>  50 

n  ■ 

>  100 

n  ■ 

>  200 

MLE 

BLUE-k 

MLE 

BLUE-k 

MLE 

BLUE-k 

MLE 

BLUE-k 

0.1 

1.53 

.52 

1.32 

1.02 

1.19 

1.10 

1.08 

1.09 

0.2 

1.46 

.95 

1.26 

1.17 

1.13 

1.24 

1.08 

1.30 

0.3 

1.35 

1.15 

1.16 

1.22 

1.08 

1.28 

1.03 

1.35 

0.4 

1.28 

1.09 

1.15 

1.24 

1.03 

1.30 

1.03 

1.43 

O.S 

1.28 

1.18 

1.05 

1.20 

1.04 

1.36 

1.03 

1.45 

0.6 

1.16 

1.10 

1.11 

1.30 

1.01 

1.38 

1.02 

1.49 

0.8 

1.29 

1.22 

1.11 

1.32 

.99 

1.36 

1.02 

1.52 

1.0 

1.26 

1.23 

1.09 

1.32 

1.06 

1.43 

1.00 

1.S4 

1.2 

1.17 

1.24 

1.08 

1.35 

1.04 

1.45 

1.00 

1.53 

1.6 

1.22 

1.31 

1.01 

1.29 

1.04 

1.52 

1.00 

1.56 

2.0 

1.08 

1.17 

1.11 

1.46 

1.06 

1.52 

.99 

1.54 

lie  Standard  Deviations 


n 

•  25 

n 

-  50 

n 

•  100 

n 

-  200 

1  „ 

MLE 

BLUE-k 

MLE 

BLUE-k 

MLE 

BLUE-k 

MLE 

BLUE-k 

_ 

1.22 

.99 

.78 

.58 

.55 

.40 

.42 

0.2 

1.55 

.68 

.80 

.81 

.44 

.55 

.30 

.50 

m 

1.06 

.86 

.51 

.63 

.35 

.52 

.21 

.48 

Iran 

.68 

.56 

.67 

.34 

.53 

.20 

.56 

HI 

.85 

.15 

.47 

.54 

.30 

.62 

.17 

.54 

m 

.74 

.81 

.44 

.58 

.27 

.58 

.18 

.58 

0.8 

.78  .66 

.47 

.64 

.21 

.52 

.16 

.61 

111 

1.18 

.90 

.42 

.65 

.28 

.56 

.17 

.63 

IQBBI 

.39 

.65 

.23  .59 

.14 

.59 

im 

.75 

.87 

.28 

.50 

.22 

.65 

.16 

.64 

2.0 

.48 

.69 

.47 

.79 

.23 

.64 

.13 

.59 

Table  6.  Simulation  experiments  (100  realizations)  comparing  estimates 
of  the  mixed  exponential  distribution  function  by  joint  MLE 
estimates  of  parameters  and  BLUE-k  estimate  of  S  with  known 
a  .  Type  I  censoring  of  data. 


Average  of  Sum  of  Absolute  Differences  from  true 
distributions . 


a 

MLE 

BLUE-k 

MLE 

BLUE-k 

MLE 

BLUE-k 

MLE 

BLUE 

0.2 

.49 

.62 

CM 

K1 

. 

.34 

.23 

.24 

.17 

.16 

0.4 

.45 

.59 

.32 

.38 

.24 

.23 

.14 

.13 

0.6 

.41 

.55 

.35 

.37 

.21 

.21 

.14 

.15 

0.8 

.41 

.46 

.27 

.27 

.22 

.23 

.14 

.14 

1.0 

.40 

.47 

.27 

.31 

.24 

.24 

.IS 

.15 

l.S 

.41 

.48 

.34 

.32 

.21 

.21 

.16 

.15 

2.0 

.39 

.54 

.33 

.34 

.22 

.22 

.13 

.13 

2.S 

.40 

.48 

.33 

.37 

.21 

.20 

.13 

.12 

3.0 

.38 

.51 

.31 

.31 

.21 

.19 

.14 

.13 

Table  7.  Simulation  experiments  (100  realizations)  comparing  estimates 
of  the  mixed  exponential  distribution  function  by  joint  MLE 
estimates  of  parameters  and  BLUE-k  estimate  of  S  with  known 
a  .  Type  I  censoring  of  data. 


0.8 


Average  of  Sum  of  Absolute  Differences  from  true 
distributions 


MLE.  BLUE-k  MLE  BLUE-k  MLE  BLUE-k  MLE  BLUE-k 


n 

•  25 

MLE 

BLUE-k 

_ 

.68 

.54 

.52 

.48 

.53 

.49 

.51 

.51 

.51 

.49 

.53 

.51 

.53 

.53 

.42 

.41 

.43 

.43 

.31  I  .22 


.32  .21 


.25 


.22 


n 

•  200 

MLE 

BLUE-k 

.21 

.18 

.18 

.16 

.15 

.15 

.14 

.13 

.14 

.12 

.14 

.14 

.15 

.14 

.14 

.13 

.15 

.15 

21 


Table  8 


Simulation  experiments  Cl 00  realizations)  comparing  estimates 
of  the  mixed  exponential  distribution  functions  by  joint  MLE 
estimates  of  the  parameters  and  method  of  moments  of  the 
parameters.  Complete  samples. 


Average  of  Sum  of  Absolute  Differences 


n  » 

25 

n  » 

50 

n  ■ 

100 

n 

.. 

«  200 

a 

MLE 

MofM 

MLE 

MofM 

MLE 

MofM 

MLE 

MofM 

2.S 

.27 

2.17 

.18 

1.45 

.14 

.76 

.10 

.52 

3.0 

.26 

1.38 

.18' 

.83 

.14 

.67 

.09 

.33 

3.5 

.27 

1.06 

.19 

.78 

.13 

.74 

.09 

.51 

4.0 

.22 

.79 

.18 

.87  |  .14 

.60 

.09 

.60 

